General Disclaimer 


One or more of the Following Statements may affect this Document 


• This document has been reproduced from the best copy furnished by the 
organizational source. It is being released in the interest of making available as 
much information as possible. 


• This document may contain data, which exceeds the sheet parameters. It was 
furnished in this condition by the organizational source and is the best copy 
available. 


• This document may contain tone-on-tone or color graphs, charts and/or pictures, 
which have been reproduced in black and white. 


• This document is paginated as submitted by the original source. 


• Portions of this document are not fully legible due to the historical nature of some 
of the material. However, it is the best reproduction available from the original 
submission. 


Produced by the NASA Center for Aerospace Information (CASI) 






















iHi5l (0) PARTICLES. ENV/I « pmmpmtc 

'« ?HE JcJ??r N^6-26ii6 

ATMLSPHERE (CORNELL UNIV.) 69 P HC »4.50 

CSCL 03B UNCLAS 

G3/91 42232 


CORNELL UNIVERSITY 

Center for Radiophysics and Space Research ^ 


ITHACA, N.Y. 


CRSR 617 


PARTICLES, ENVIRONMENTS AND POSSIBLE ECOLOGIES 

JOVIAN ATMOSPHERE 
(Revised Version) 


Carl Sagan and E. E. Salpeter 


(s 


'^L9 5li 







PARTICLES, ENVIRONMENTS AND POSSIBLE ECOLOGIES 
IN THE JOVIAN ATMOSPHERE 

I 

Carl Sagan and E. E. Salpeter 

i 

Original received: November 1975 

Revised version received: May 1976 


Center for Radiophysics find Space Research 
Cornell Dnivert’ity 
Ithaca, N.Y. 14853 


r 

i 

I 














Abstract 


The eddy diffusion coefficient is estimated as a function of alti- 
tude, separately for the Jovian troposphere and mesosphere. The 
growth-rate and motion of particles is estimated for various substances: 
the water clouds are probably nucleated by NH^Cl, and sodium compounds are 
likely to be absent at and above the levels of the water clouds. Complex organic molecules 
produced by the Ly a photolysis of methane may possibly be the absorbers in the 
lower mesosphere which account for the low reflectivity of Jupiter in 
the near ultraviolet. The optical frequency chromophores are localized 
at or just below the Jovian tropopause. Candidate chromophore molecules 
must satisfy the condition that they are produced sufficiently rapidly 
that convective pyrolysis maintains the observed chromophore optical depth. 
Organic molecules and polymeric sulfur produced through H 2 S photolysis at A > 

2300 Aprobably fail this test, even if a slow, deep circulation pattern, driven 
by latent heat, is present. The condition may be satisfied if complex organic chromo- 

O 

phores are produced with high quantum yield by NH 3 photolysis at X < 2300 k. However, 
Jovian photoautotrophs in the upper troposphere satisfy this condition well, even with 
fast circulation, assuming only biochemical properties of comparable terrestrial 
organisms. Unless buoyancy can be achieved, a hypothetical organism 
drifts downwards and is pyrolyzed. An organism in the form of a thin, 
gas-filled balloon can grow fast enough to replicate if Ci) it can sur- 
vive at the low mesospheric temperatures, or if (il) photosynthesis 
occurs in the troposphere. If hypothetical organisms are capable of 
slow, powered locomotion and coalescence, they can grow large enough 
to achieve buoyancy. Ecological niches for sinker's, floaters and 
hunters appear to exist in the Jovian atmosphere. 


I. INTRODUCTION 


The contemporary Jovian atmosphere has some similarities to the 
primitive terrestrial atmosphere, where organic molecules were pro- 
duced readily; and experiments in which a wide variety of organic 
molecules are produced under conditions similar to those of the Jovian 
atmosphere have been carried out (Sagan and Miller, 1960; Sagan et al. , 
1967; Woeller and Poniiamperuma , 1969; Rabinowitz et al, , 1969; Sagan 
and Khare, 1971a, 1971b; Khare and Sagan, 1973; Ferris and Chen, 1975). This has 
led to the hypothesis that there may be a Jovian biology (Sagan, 1961; 

Shklovskii and Sagan, 1966). Observations in the visible and in the 
ultraviolet have implied the presence of particulate matter in the 
Jovian atmosphere, which may (but need not necessarily) be connected with 
questions of Jovian organic chemistry and biology. Both chemical and biological issues 
are affected by fluid dynamics in the Jovian atmosphere. We are here 
concerned largely with these effects. 

Modern models of the structure of Jupiter (for references, see, 
e.g., Stevenson and Salpeter, 1976) predict that the bulk constituent, 
a hydrogen-helium mixture, is fluid throughout with temperatures in- 
creasing inx-7ards to above 10^°K. Plausible organisms require tempera- 

3 

tures well below 10 One principle problem which must be faced by 

such hypothetical organisms is descent to pyrolytic depths due either 
to the acceleration of gravity or to convective downdrafts. The second 
half of the present paper addresses (1) whether such organisms unde' go 


~z~ 


significant growth and replication in less than the timescale for drift " ; 

to pyrolytic depths; or (2) whether other adaptations to avoid pyrolysis 
are hydrodynamically and hiologically feasible. We shall see that 
(1) becomes easier for smaller organisms and (2) for larger ones. 

That plausible and internally consistent Jovian ecologies can be 
described cannot , by itself, demonstrate the likelihood of life on 
Jupiter. A tenable argument for life on Jupiter must also show that ; ! 

the origin of life on that planet was possible. Since many steps, 
particularly the later ones, which led to the origin of life on Earth 
are only incompletely understood (see, e.g., Miller and Orgel, 1974) a 
thorough discussion of the origin of life on Jupiter cannot be under- 
taken at the present time. Host skepticism on this subject seems to 
arise from the contention that synthesized molecules will be carried 

* J 

convectively to pyrolytic depths before the origin of life can occur, | 

The unspoken premise in such an argument is that pyrolysis would not ■ j 

• 1 

have occurred on the primitive Earth. However, the measured Arrhenius | 

rate constants imply through a regression analysis that many essential 

■'Vi 

organic molecules, for example a number of biological amino acids, suffer ] 

thermal degradation at typical terrestrial surface temperatures in ' i 

■ ! i 

\ 

geologically short time periods. The very limited available experi- . .j 

o ' 

mental data show^ for example, that the half life at 30 G of the 

4 5 ^ 

amino acid serine is y 10 yrs and of threonine is < 10 yrs; while at 

the same temperature the half life of the simplest amino acid alanine 

is >10^^ yrs (Vallentyne, 1964). Because of the exponential Arrhenius 

kinetics, the highest temperatures on a temperaturer-variahle planet 

dominate the net molecular degradations; also there is evidence 

i 

that aqueous solution and alkaline pH — - both expected on the 
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primitive Earth — enhance the decomposition rates (Vallentyne, 

196A). It does not seem at all improbahle that some essential 

precursor molecules for the origin of life had half lives on the 

2 3 

primitive Earth of 10 to 10 yrs, not greatly dissimilar to the 
lifetimes against convective pyrolysis of organic molecules on 
contemporary Jupiter, as calculated below. 

Many workers on the origin of life have suggested that the 

rate-limiting steps were stochastic in nature. The volume available 

■ ■ 3 

for natural molecular experiments on Jupiter is approximately 10 
times larger than on Earth. If phase interfaces play a critical 
role in the origin of life (see, e.g., Bernal, 1973) Jupiter — with 
a great concentration of solid and liquid atmospheric aerosol 
particles — - would be particularly favorable. 

The origin of life on Earth took a period < 0.1 the lifetime 
of the Earth, and possibly much less (Sagan, 1975). On Jupiter, 
we have much longer periods of time available for chemical experi- 
ments, a larger interaction volume, and a much larger area for 
phase interface chemistry. On the other hand the thermal degrada- 
tion timescales are shorter for many molecules, and there is a much 
greater excess of molecular hydrogen than on the primitive Earth, 
where H 2 was rapidly depleted during planetary formation and in 
subsequent Jeans escape. The effect of molecular hydrogen on pre— 
biological organic chemistry is inadequately understood, but is briefly 
discussed below. Its effect on quenching gas phase organic reactions 
is less than linear, at least for some reactions. We conclude that the absolute 
a priori likelihood of life on Jupiter, as on Earth, cannot be 
estimated reliably at the present time, but * — as on Earth— - can 
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hardly be ruled out. The relative a priori probabilities of the origin of life, in 
a comparison of Jupiter and Earth, depend on which reactions and molecular precursors 
are rate-limiting.. If, for example, the alanine abundance is rate-limiting, the 
Earth is strongly favored; if serine is rate-limiting, Jupiter may be favored. It 
is also possible that there exist a variety of different paths to the origin of life, 
for Si;me of which the Earth is favored and for others of which Jupiter is favored. 

We devote the remainder of the present paper to a characteriza- 
tion of the relevant physical environment of Jupiter; to a discussion 
of the chroraophores responsible for the observed coloration of the 
planet; and to an analysis of some permissible ecological niches of 
hypothetical Jovian organisms. 

Values for the eddy diffusion coefficient K play an important 
role in such arguments* Estimates are given in Section II, separately 
for the convective troposphere and for the stabler mesosphere. 

Analytic approximations are also given there for abundance as a function 
of height for particles which are produced at some upper level and are 
transported downwards towards pyrolysis. 

We are concerned with two sorts of particulate matter: particles 
that condense out under conditions of thermal equilibrium and particles 
that require an investment of free energy for their production. 
Equilibrium condensation has been studied previously for H2O, WH^SH 
and we investigate it further in Section HI for a number of 

substances with higher condensation temperatures. We estimate 
growth-rates for these condensates and discuss whether NaCl or other 
particles can be present at the level of the water clouds to act as 
nuclqation agents. The production of chromophores and of non-equilibrium 
organic molecules is treated in Section IV. Wd obtain a lower limit to 
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the production rate of ethane from methane photolysis by solat 
Lyraan-alpha; and treat organic chromophores produced by various 
wavelengths of light, as well as polymeric sulfur chromophores. 

The motion of hypothetical organisms is discussed in Section V; 
the diffusion of metabolities and the consequent growth of organisms 
in Section VI. We consider four different cases depending on 
whether organisms can or cannot survive at the low temperatures of 
the mesosphere and whether they can or cannot utilize visible light 
for photosynthesis. In Section VII we consider the more rapid growth 
attainable if organisms can control their motions. 


II. EDDY DIFFUSION IN THE JOVIAN ATMOSPHERE 

The total infrared flux emerging from Jupiter (Aumann, et al. » 1969; 
Ingeraoll, at al , , 1975) is = a ^eff 

^"^ef f ~ 127*^K) , with about half of this flux coming from the deep 

4 -2 -1 

interior ~ 0-8 x 10 erg cm s ). The radiative opacity in- 

creases fairly rapidly with increasing pressure, P, and radiative-con- 
vecti »e model calculations predict a fairly rapid change across the 
boundary level defined here as the tropopause. Above the tropopause 
we have radiative C subadiabatic) heat transport in the mesosphere; 
below we have convective heat transport (with the temperature gradient 
only slightly superadiabatic) in the troposphere. Model atmospheres 
(Trafton and Munch, 1969; Hogan, et al . , 1969; Wallace, et al . , 1974) 
and spectroscopic observations (Gillett, et al . , 1969; Aitken and 
Jones, 1972; Lacy, et al . , 1975) suggest a temperature T ~ 175*^K. at 
the pressure level of P = 1 bar and place the tropopause somewhere 
slightly above this. We shall adopt the level (T 145“, P 0.5 bar) 
for the tropopause. These numbers are consistent with Pioneer 10 
infrared radiometer data (Orton, 1975a) and with ground-based Infrared 
radiometric measures (Orton 1975b) and possibly also with Pioneer 10 
and 11 S-band occultation data (Kliore and Woiceshyn 1976). We assume 
a helium abundance of about 10 percent by number, consistent with separate 
measurements from the occultation of Beta Scorpii (Elliot, et al. , 

1974) ; from Pioneer 10 and 11 ulcraviolet photometry (Carlson and 


6 



7 


Judge, 1974); and from the pressure hroadening of the J = 0 and J = 1 
hydrogen lines (Houclc, et al. , 1975.) Assuming a helium-abundance of 
about 10% (by number), it then follows that y/CY - 1) * 3» where y is 
the ratio of specific heats, and that the mean molecular weight y " 2.2. 
For pressure P, pressure scale height H = P (dP/dz) , density p and 
kinematic viscosity v in the troposphere we adopt 


P « 

p « 
V a 


^175®K^ ® ™ ^1750K^ 27km; 


T 2 —4 —3 

(•JT^) 1.6 X 10 gm cm , 


( 


— 0.5 cm^ s”^ 
175°K' 


( 1 ) 


In the mesosphere, the temperature gradient decreases rapidly (for 
P < 0.2 bar, say), and the minimum temperature is about 120°K. fol- 

lowed by a temperature inversion whose structure is not yet well-known 

(possibly T ** 150°K at P n. 10 ^ bar) . 

^ max 


For a general convective, circulation pattern the vertical con- 
vective heat-flux F^ is given by 

^c f2) 

where c is the specific heat at constant pressure, AT is the typical 
vertical excess in the temperature drop and w^ the typical vertical 
convection speed. For vertical distances larger than the Vertical 
extent A of the circulation pattern one has an eddy diffusion coef- 
ficient w^A and we shall denote hyt the eddy diffusion time, 

H/K^, for diffusion through a single pressure scale height H. In 
the simplest form of Prandtl mixing-length theory (see, e.g., 

Schwarzschild, 1958) the circulation is characterized by two parameters 

O' 

the mixing-length A and the superadiabaticity with w^ aiid AT/T both 
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proportional to the superadiabaticity. We replace F(, by 1.4 x 10 erg cm s , a 
typical value for the Jovian troposphere, employ eq. (1) and write K . w and t 

in the form 

= C (T/175“K)^^^ X 10^ cm^ s“^ , 

Wc “ C^^^(T/175°K)“2/3 X 76 cm s“^ , (3) 

t^ = (T/175'’K)^^^ X 7.5 X 10^ s . 

In the notation of mixing-length theory our parameter 5 equals 

A / 3 

(A/0.05H) ^ , The Jovian tropospheric lapse rate is very close 

to adiabatic (AT/T 10 ), as is true by the same argument (Sagan, 

I960) for all other planetary tropospheres. 

In the simplest version of mixing-length theory H and ^ 100. 

The actual circulation patterns in the troposphere are likely to be 

much more complex (Hunten 1975). In particular one should consider 

the effects of rotation at great depths and of latent heat release in the 

cloud layers near ^:he tropopause (Gierasch, 1973), For a constant value 

of the heat-flux F^ in eq. (2), these effects usually increase AT and 

decrease w considerably, leading to smaller values of K and ? in 
c c 

eq. (3). These considerations suggest , as a simple working hypothesis, ^ = 1 which 
should be accurate to within a factor of about 200. 

Gierasch (1976) has suggested a specific large-scale circulation pattern which 
extends all the way from the bottom of the water-cloud layer (T 'u 280®K) to the tropo- 
pause. This pattern is driven by the latent heat of water which leads to a relatively 

-2 

large value of AT/T 10 . This corresponds to ^ O.pl with 

timescales for complete turnover of the region 
of up to a few years, even though some cloud patterns ( above the water 

clouds) are observed to change in shorter times (Peek, 1958). On 

Gierasch* s model, the eddy diffusion coefficient increases by a very . 

large factor just below the cloud-bottom to ? 'u 10. 
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In the mesosphere the bulk of the host flux is carried by 

radiation, but some fraction of the flux could drive a circulation 
pattern working against a subadiabatic temperature gradient. The 
vertical speed of such a circulation pattern must satisfy an in- 
equality, based on the fact that the adiabatic cooling rate cannot 
exceed the total available heating rate: 

Hpw^ [Cp(dT/dz) + g] < (4) 

Except in the immediate vicinity of the tropopause, the temperature 
gradient in eq . (4) can be neglected for the troposphere* Putting 
T - 150°K and assuming that the vertical length-scale for the circula- 
tion pattern is of order H, the eddy diffusion coefficient K is of 

m 

order w H, so that 
m 

w < (1 bar/p) 0.016 cm s K < (1 bar/P) ^ x lO^ cm^ s (5) 
in — m. — 

Thus, typical velocities just below the tropopause are many thousands 
of timeslarger than just above the tropopause. Analysis of the ob- 
served by a albedo of the upper mesosphere (Wallace and Hunten, 1973, 

Carlson and Judge, 1974) indicates K 3 x 10 ~ cm s at levels 

m 

-7 

where P 3 x 10 bar; analysis of B Sco occultation measurements 

6 2—1 

(Sagan et al. , 1974) indicates > 10 cm s at similar levels. 

It is theoretically plausible (French and Gierasch, 1974) that gravity 
vrava propagation upwards in fhe mesosphere con f r I'hut ea lu the eddy 

tllfi'usiou. t'or a single wave mode, has a step-f uiic L Jon behavior, 
but overall it seems plausible to assume .;ome constant fraction of 
the maximum expression in eq* (5). The value given by the Ly ® albed 

then suggests 

(1 bar/P) 100 cm^ (6) 
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uncertain by factors of order 10 or 100. Note that K Tem- 

m 

perature T and K(z) are plotted against pressure, schematically, in 
Figure 1 , [The dashed curve for K corresponds to Gierasch’s (1976) 
model in which the latent heat of water condensation is taken into 
account.] 

We shall need continuity equations for n^(z), the number density 
as a function of height z of some complex molecule (of molecular 
weight >> 2.2) or of some particle. Let w^(z) be the steady-state 

vertical drift velocity of the molecule or droplet falling under 
gravity through an atmosphere with pressure scale height H, eddy 
diffusion coefficient K and mean molecular weight 2.2. In standard 
notation, the continuity equation (see, e.g,, Strobel, 1973) reads 

(K + D ) A + .SL) n = -f , 

'' dP^ ^ ^ °j j ’ 

where is the effective molecular (or particle) diffusion coefficient 
and H. = H (2.2/p.) << H is the effective scale height for the j th 

j 1 

particles alone. ®. is the vertical flux of particles of species j. 

3 

and the drift-velocity for free fall under gravity are related 

by D. = w,H.. We shall be interested only in levels well below the 
J 3 

turbopause where D. « K and give only a simplified version of the standard treatment 
3 

(Strobel 1973; Prinn 1973) . The steady’-state continuity equation then reduces to 


K(z) + w^Cz)] n^(z) - -^.Cz) 


In Section IV we shall apply eq . (7) to the number density 

n(z) of organic molecules or other chromophore particles for which 


the drift-velocity can be neglected (w^ << K/H) • Assume further that 
the particles are mainly produced at some higher level z 13 > 0 

and destroyed at some lower level z = Q, so that the flux $ equals 
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the total production rate (cm s ) for 0 < z < and » 0 outside 

that range. If the production level is in the troposphere, the 

eddy diffusion coefficient is given fay in eq . (3), which varies 

quite slowly over a scale-height H. An approximate solution 

of eq. (7), which neglects terms in dK /dz, is then 

c ' 

n(z) S (§H/K ) (1 - for 0 < z < Z , 

o 

« (OH/K ) [ “ ^0^^® - e"^o^®] for Z < z. 

Thus, particle number densities decline exponentially with altitude 

above the production level and decrease more slowly with depth below 

the production level. Mathematically, the solution in eq. (8) depends 

on the boundary condition that n(z) = 0 for z ^ 0 and that $ is already 

constant for infinitesimally small positive z. Physically, the require 

ment is merely that the destruction rate change from very small values 

(at positive z) to very large values (at negative z) over less than a 

scale height. This should be a good approximation since pyrolysis 

rates are very temperature sensitive. (But note that the level z = 0 

refers to different altitudes for different molecular species.) 

The situation is more complicated if the production level Z^ is 

in the mesosphere and the destruction level (z H 0) in the troposphere, 

with z = Zp < Z^ representing the tropopause. However, neglecting 

again derivatives of K = K (for z < Z ) and assuming that K = K 

c p “ m 

25 /TT 

(for Z^ < z) is proportional to e as in eq, (6), the continuity 

eq. (7) can again be solved analytically. With the fact that 

Z >> E and K (Z ) >> K (Z ) for cases of practical interest, the 
p c P m p 

solution can be approximated by 
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. I 

I 
! 

i 

f 
i 

•J 
\ 

I n(z) w (§E/K^) (1 - for 0 < z < Z^, 

i ■ ' 

M > Z + H £or Z < z < Z^, 

i HI c p _ 

I ® for < z. 

■ I , 

'•I 

Thus, particle number densities increase slowly with 
altitude above the destruction level until the tropopause; decline 
slower than exponentially (“ z e between the tropopause and the 

production level; and decline exponentially above the production level. 
For mesospheric production and tropospheric destruction , we therefore 
I expect maximam rtuiaber densities just above the tropopause* Againj no 

j particles exist at Z < 0. Note that the mixing-ratio (“^^/pgas) remains 

j constant above the production level whether this level is in the 
f mesosphere or troposphere. In reality, this mixing ratio would decrease 

[ slightly with increasing height z if some additional destruction takes 

I place above the production level or if the drift-velocity is not 

I . .. ■ 

i entirely negligible. 

■■ i . . - - ^ ‘ - 

. i . . ■ ■ ■ 

i ■ . . 

f ' ' ■ 

f . 




III. THERMAL CONDENSATION IN DEER CLOUD LAYERS 

The thermal condensation of droplets to form cloud layers in 
the Jovian atmosphere has heen discussed (Weidenschilling and Lewis, 
1973; Stauffer and Kiang, 1974) only for the upper levels of the 
troposphere (T S 300°K) , where direct observations are possible. 
Molecules involving H, C, N', 0, P and S have been considered, which 
lead especially to condensation of H 2 O (T ~ 275°K) , NH^SH <T ~ 220^K) 
and NH^ (T ~ 160*^K) . The condensation of droplets of more refractory 
materials in deeper layers (higher temperature, T) is also of some in- 
terest, especially for compounds containing Na and/or Cl. If such 
small grains or droplets can be convected to the upper troposphere, 
they can act as seed-nuclei for the condensation of H 20 -drops 
(at T 5 ; 280°K) . Their abundance is also relevant for any hypothetical 
Jovian biology and for the question of the replication of terrestrial 
microbial contaminants in the Jovian clouds. The following crude 
estimates and conjectures should serve at least for purposes of 
orientation. 

In the very deep interior of Jupiter where the density exceeds 
0.2 gm cm ^ (P > 10^ bar, T ~ 10^°K) , incomplete solubility of var- 
ious substances in hydrogen Could lead to substantial chemical 
separation. This possibility, discussed so far only for He (Salpeter, 
1973; Smpluchowski, 1973), should be investigated further; however, 
we assume provisionally cosmic abundances of all elements (Allen, 1973) 

■ . 13 , : . 
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at P 10^*^ bar (T 2,500°K), ,say. At this level practically all 
substances are gaseous. Refractory molecules begin to condense at 
slightly lower temperatures. 

For each material, with £ the cosmic abundance by number (mixing 
ratio) of the least abundant element in the compound, we must evaluate 
a number of quantities; 

(1) The equilibrium condenaatlGn temperature appropriate 

to the Jovian troposphere, i.e., that temperature at which 
the saturation vapor pressure equals f times the pressure 
in eq. (1). Values for have already been given by Lewis (1969b). 

(2) Assuming the absence of any external seed-nuclei, we cal- 
culate the amount of supercooling required from 

homogeneous nucleatipn theory. For this purpose w& need a parameter 
ri, defined in eq. (7) of Salpeter (1974), which is essentially the 
number of times appropriate molecules of a particular species stick to a 
given surface site on a droplet during a typical flow-time. In that 
definition for q we use s A/ 0,5 for the surface sticking coefficient, eq. (1) 
for the pressure and replace the flow-velocity v by the typical 

convection speed x<r^ in eq. (3), which gives q 'v f (T^^^/175'’K)^*^ 

■ ^ ^ 13 ■ 

X 3 X 10 . Values of f, and T (taken from JANAF 1971) are 

con 

' ' ■ ’ 8 ■ ■ 

given in Table 1 for a number of compounds and T] exceeds 4 x 10 

in ail cases. The ratio AT /T is given by eq. (9) of Salpeter 
(1974), but multiplied by a factor 6 ^ (correcting an error kindly 
pointed out by B. Draine) , This ratio is fhirly small in the 
cases considered and not of direct importance for our crude estimates. 

But we do need the final radius a^ of typical drops which would 
result at the end of purely homogeneous nucleation. 
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At first sight one would expect to be of order n times the 

intermolecular spacing in the liquid, but there are a number of 

correction factors. The most important correction is due to the fact 

that the drop-radius a becomes large compared with the mean free path 

—2 —5 

for a molecule in the gas, (T/175“K) x 6 x 10 cm. The 

droplet's growth is controlled by the diffusion speed rather than 
the gas-kinetic speed and the rate is multiplied by a Correction 
factor (JZ-^/a) « 1. Including some further logarithmic factors we 
find, as a very rough approximation, 

a^/u (i75“K/T)n^^^(log3_QTi)"^'^^ X 10"^ cm . (10) 


Our values for a^ are given in Table 1 and have to be compared with 
a number of other length scales; 

(3) Let a be a critical radius for a spherical droplet such that the 

hydrodynamic drag-coefficient is unity. For a « a^^^ the Reynolds 
number is small and the drag-coefficient is given by the Stokes for- 
mula Gjj = (12v/aw^), where is the Stokes expression for the terminal 
vertical drift-velocity of the sphere. Using eq. (1) we have 



0.007 cm 



175°k | 

T J 



Wd(a) 


2 

9 


2 

V 


dr 


a. 


0.15 cm s 


-I 


l-% 


il75“K 



_3 

vjhere is the internal density (in gm cm ) of a liquid drop, 
pother critical length is 0^, the droplet radius for which the drift- 
velocity w^ equals the typical convection speed in eq, (3). For 
the cases of interest (T < 2,500°K) we have ad*' so that eq. (11) 
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holds. With ^ 2 and ? 1 in eq. (3) we have 

0.002 cm (T/175“K)"^^^^ (12) 

Large drops with a » a^ ’’rain out" rapidly, whereas small droplets with a « 
are carried along by the convection pattern and can move up to higher levels. 

For all cases in Table i (except possibly for KH^Cl) we have a^ > a^, so 
that homogeneous nucleation leads to large enough drops for most of the matter 
to be rained out. However, if a very large number of nucleation seeds are 
present originally, a large number of much smaller droplets condense out and 
the droplets must grow before they can fall out as rain. Let that radius 

for which Cjj = *^dr^^g* 

^coa X 3 X lO^^cm (13) 


Once a ^ a , droplets usually can grow very quickly by coalescence, 

but for a < a growth must occur by the slower process of smaller 
coa 

drops evaporating and the excess vapor diffusing to larger drops 

(Langmuir, 1948; Mason, 1972; Rossow, 1976). An important parameter 

n is the growth rate from this process (from radius 0,5 a to 

vap coa 

2 

a , say) expressed in units of H /K , a typical convection time, 
coa . C- 

With an uncertainty of at least a factor of 10, we have 


H '\< 4 X 10 f (- 
vap 


175 


7/6 


x^y' 


(14) 


If n » 1 for a particular material, most of it will form 

■ ■ V sp 

into large "rain-drops" in a cloud layer at levels close to the. con- 
densation temperature, and only small amounts of the material 

will be found at higher, cooler levels. If t) << 1 

vap 
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one might expect only small droplets of the material to form, so that 

w, « w and upward convection can maintain an appreciable mixing 
d c 

ratio of the droplets at higher levels. In practice, however, 
various complications usually tend towards more 

-1 ' 

rapid formation of larger droplets (Rossows 1976). For 10 

or 10 the circumstances are probably mixed: Some rain near the 

T level, some small droplets a few scale heights higher up, but 
con 

with a mixing ratio which decreases slightly with increasing height. 

Values for from eq. (14) are given in Table 1. Except possibly for NH^^Cl, 

-2 

n » 10 as well as a_ > aj and most of the material rains out and cannot 
vap f d 

move up to higher levels. However, in upward moving convection cells the 

saturation vapor pressure continues to fall; droplets then continue to condense 

out, but with decreasing values of the mixing ratio f and hence of the parameter 

f|. Let f^ be that (decreased) value of the mixing ratio, for which the droplet 

size a^ (due to homogeneous nucleation) in eq. (10) equals a^» Values of are 

also given in Table 1. The mixing ratio will not drop appreciably below f^ in 

upward moving convection cells, since droplets are now too small to rain out 

rapidly enough (and further groirth of droplets is also too slow). 

For the cases with T 5- 1600'’K the radii satisfy a- > a , coalescence 

con . f coa 

is rapid and large drops are rained out immediately. At the high pressures in 
Jupiter's interior, liquid Mg^SiO^ (as well as Ca- and Al~ silicates and metallic Fe) 
Condenses before MgSiO^ and before Si02* When liquid Si02 begins 
condensing near the 1800°K-levei, most of the metallic elements 
have already been rained out, but sodium has not. For pure Ha 2 SiO^ 
our estimates indicate a slightly lower condensation temperature than 
for Si02> which introduces an uncertainty; It is conceivable that all 
the silicon, has been rained out before sodium condenses, in which case 
gaseoUs NaCJZ- and HaOH (since the elemental abundance of Na exceeds 
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that of CA) survive to higher, cooler levels. It is more likely that 
much of the sodium condenses out between 1800°K and 1600^ in the 
form of some sodium silicate-silica solution, in which case gaseous 
HCA and some NaC^ survive to higher levels. In either case, more of 
the surviving NaCit condenses out near 1000°K; in the more likely case, 
KH^CA condenses near 400°K; in the less likely case Na 2 S condenses 
near lOOO^^K. 

We have some tentative implications for levels with T ^ 700°K, 
where Jovian biology (or at least organic chemistry) might be possible; 
Sodium (or other alkali metals and alkaline earths) and chlorine (or 
other halogens) cannot be present simultaneously. However, either 
sodium sulfide or (more likely) ammonium chloride is present. For 
the more likely case, some fraction of the NH^C£ - crystals are con- 
vected upward and act as seed-nuclei for the condensation of water 
droplets, so that essentially no supercooling is present in the base 
level of the Water clouds (Stauffer and Kiang, 1974) • For the less 
likely case of Na 2 S crystals convection over a larger height dif- 
ference is required for nucleating the water, and is slightly 
smaller. However, the chromophores to be discussed in Section IV, 
as well as Na 2 S and NH^Cl, could probably also act as seed-nuclei, 
and it seems generally safe to assume no supersaturation for water. 

We estimate below the amount of exogenous sodium which may arrive 
at the cloud level from micrometeoritic inf all. 

Since both sulfur and sodium ha"ve been detected in toroidal nebulae 
about Jupiter - — probably arising from the surface of lo — these 
clouds are a potential source of S and N a in the clouds of Jupiter. 
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Ftom observations of forbidden SII circumjovian line emission, Brown 
(1976) calculates a total sulfur flux from the surface of lo of 
10 cm sec , squxvalent to roughly 10 cm sec at the 

Jovian clouds. This is some eight orders of magnitude less than that 
expected from H^S photodissociation (see Section IV), and therefore 
is entirely negligible. The cosmic abundance of Na is at least an 
order of magnitude less than that of S, and so only very small quantities 
of sodium could be present at the Jovian clouds due to leakage from 
the lo Na toroid. 

Approximate influxes of micrometeoritic sodium and sulfur, 
assuming the fluxes are the same on Jupiter as on Earth and that 


micrometeorites are primarily cometary ices with cosmic abundances 


>5 -2 -1 


6 -2 -1 


of other atoms, are *''10 cm sec for Ka, and cm ** ge 


c tor 


S. Thus, the micrometeoritic infall rate of these atoms is likely to 
exceed that from the lo-associated toroidal clouds . In the upper 


troposphere with T ^ 175®K, eqs. (3) and (9) then give n 10 Na 

-3 -15 

atoms cm , corresponding to a mixing ratio f ^ 10 , and even 

smaller at the water cloud levels. These values are small compared 

with the values of f in Table 1. 

m 
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TABLE 1 

The equilibrimn condensation temperature, and estimates of the relevant 

mixing ratio, f, the radius a_, the parameter T) and the minimum mixing ratio 

t vap 

f for a Tiumber of substances above their condensation levels* 
m 



Fe 

MgjSiQ^ 

SiOj 


: =3: 

NaCl 

NH4C1 

(Na^S) 

-IbgioE 

4.1 

4.6 

4.4 

5.7 

6.3 

6.3 

(5.7) 

T /lOQ^K 
con 

25 

23 

18 

16 

10 

4 

10 

a--/ 10 ^ cm 
t 

120 

60 

40 

20 

4 

2 

8 

r\ 

vap 

60 

20 

2 

1 

0.1 

0.03 

0.4 

-lOgn 

^10 m 

7.9 

7.7 

7.2 

7.8 

7.0 

6.4 

7.0 


IV, PHOTOPRODUCTION OF MOLECULES AND CHROMOPHORES 


(a) Far UV : 

At thermal equilibrium methane is by far the most abundant carbon- 

compound in the Jovian atmosphere, where it remains in gaseous form. Methane 

o 

absorbs UV photons at all wavelengths shorter than 1600A and the total 

flux of such photons at the top of the Jovian atmosphere (mainly from 

o 9 -2-1 

solar Ly a at X= 1216A) is ^ 5 x 10 photons cm s (global diurnal 

average) . Very little of the Ly a photons are scattered by neutral 

hydrogen (Carlson and Judge 1974) and about one-third of the photons process 

CH^ into other hydrocarbons (Strobel 1973a, 1975). The dominant product is 

9 -2-1 

ethane with a predicted production rate of R 1 x 10 ^2^6 ® 

-14 -2-1 

5 X 10 gm cm s , The ethane produced is rather stable and is destroyed 
mainly by eddy diffusion into the troposphere, followed by pyrolysis in deeper 
hotter layers. For an assumed mesosphere eddy diffusion coefficient K^(z) , 
accurate numerical calculations for the number density n(z) of OgHg can 
be carried out (Strobel 1975). However, eqn. (9) is sufficient for pur- 
poses of orientation (the largest uncertainty comes from uncertainties in 

B^) . Assuming the tropopause level Z to occur at P’Vp.5bar (with P 10 ^ 

—7 

to 10 bar at the production level Z ), and with K (2) given by eqn. (6) , 

o m 

the maximum ethane number density occurs one scale height H above the 
tropopause: 

n ®H/K(Z + H) /o 5 X IC?"^ cm ^ 5 X 10 ^n_ (15) 

p "^2 

where K(Z H) is the eddy diffusion coefficient evaluated at altitude 

p \ . 

Z + H. 

P 
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The maximum value of the mixing ratio occurs just below the production level 
and the total predicted column density N of ethane molecules is 

N $H^/K^(Zp) 'V' 5 X 10^^ cm~^ 3 X lO”^ gm cm”^ . (16) 

Ethane (and smaller amounts of acetylene) have indeed been detected 
in the Jovian mesosphere by means of the emission spectrum at X lOp (Ridgway 
1974) . High-resolution spectra are not yet available and the column density 
N of Q 2 Hg from the present spectra depends strongly on the distribution of 
gas temperatures, which is not well-determined by low-resolution spectra; 

If the temperature in the mesospheric inversion layer turns out to have a 

_2 

maximum value of 150°K, the observations indicate (Ridgway 1974) 0,02 gm cm 

of ethane in this high temperature region. This would be in very serious 
conflict with the predictions, since only a small fraction of the theoretical 
column density in eqn. (16) refers to the hotter region many scale heights 
above the tropopause. Presumably one is dealing ^d.th some combination of 
(a) higher mesospheric temperatures; (b) smaller diffusion rates than in 

eqn. (6); or (c) some additional source of ethane raising the production 

9 -2 -1 
rate R above 10 cm sec . 

Production rates from the same Ly a radiation in the mesosphere for 
larger, more complex organic raol6;cules has been estimated (Strobel, 1975) 
to be factors of A. 10 or 100 (by number) less than ethane, ah estimate 
consistent with the data presented by Noyes and Leighton (1941) who remark 
that the ultraviolet irradiation of ethane leads to **a surprising variety 
of products," including unsaturated polymeric hydrocarbons. Assuming 
one-tenth the rate (by mass) , the mass-production rate R is 
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5 X 10 gm cm s and eqn. (16) would give a column density for 

- 4-2 

complex organics of 3 x 10 gm cm . 

Complex organic molecules may or may not be the mesospheric chromophores 
which have been inferred from the UV-albedos: The measured albedo of Jupiter 

(Wallace et al., 1972; Savage and Caldwell, 1974) to solar photons in the 


near UV is quite low with a minimum 0.25 near X 'v 2800 A. Unit optical 

o 

depth for pure (conservative) Rayleigh scattering by H 2 for 2800 A photons 

is reached at a level with P 0.3 bar in the lower mesosphere. No simple 

molecules present near or above this level absorb in this wavelength region 

< o 

(NHj is effective only for X "v* 2,300 A) a4d complex molecules or very small 
absorbing dust grains have been invoked to provide the absorption (Axel 


1972) . Absorption in this region can be provided by a range of not im- 
plausible organic molecules (Sagan, 1968; Khare and Sagan, 1973). Amass 

° -13 2-1 

absorption coefficient near 2800 A of (2N^^g) 'v 10 cm gm would be 
required if organic molecules (beyond ethane) produced by ty a photons 
are responsible. This is within the range of measured absorption coefficients 
for complex organic molecules produced by ultraviolet light in crude Jupiter 
simulation experiments (see, e.g. , Khare and Sagan, 1973). 


(b) Cloud models 

The total thermal infrared flux emerging from Jupiter gives 

T 127"K. Most of this flux emerges at wavelengths X /v 20 ym to 60 pm; 

exf — 

the opacity sourre is mainlv il, (v/lili rtiiiftll mmiU find fv-u 
NH^); and most of the radiation emanates from levels with P "V 0.7 bar to 
0.2 bar (T 150®K to 115°K) . Only a small fraction of the thermal flux 
emerges at wavelengths X < 15 ym, but the ■ pacities of the main gaseous 
constituents of the Jovian atmosphere are low at some of these shorter 
infrared wavelengths. In the absence of any particulate matter (or complex 
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molecules), radiation with X 8»2 ym to 9.5 ym would emerge from levels 
with T 155 "K, and radiation with X 5 ym from even deeper levels with 
temperatures up to T 300®K. 

Assuming solar abundances of all elements in the Jovian interior 

and assuming thermochemical equilibrium, one can calculate condensation 

temperatures various species (cf. Lewis 1969a; Weidenschilling 

and Lewis 1973: Prinn and Owen 1976). T determines 

’ con 

the level of each "cloud bottom," and above this level the ga?^ -phase 
abundance is given by the saturation vapor pressure. The abundance and 
size-distribution of droplets (or crystals) in the upper layers of a 
cloud cannot be calculated vjithout detailed knowledge of the dynamics 
(see Sect. III). The cloud bottoms for ice (and aqueous ammonia) are 
at levels with T 280"K, for NH^SH at T 200“K, and for at T 2l 
150°K (P ^ 0.7 bar). 

Measurements of spectral reflectivities for solar radiation in the 
visible range are now available separareJ.y for the bright zones and for 
the darker belts in the Jovian equatorial and temperate latitudes (Orton 
1975c) . The fraction of the sunlight that is absorbed varies from only 
0.2 in zones to ^ 0.3 in belts for red light, but varies from 2; 0.3 in 
zones to ^0.5 in belts for blue light. The data at X 20y and 45y, 
separately for zones and belts (Orton 1975c) , are compatible with a model 
which includes an optically opaque cloudtop at the level (T ^ 145“K, 

P ^ 0i6 bar) for zones, but not for belts. This model is not unique, 
but the postulated cloudtop occurs close to the level at which solid 
crystals are expected to precipitate out if the zones represent the 
upward draft in the overall circulation. The spectral data for X ^ 8 to 14 ym 
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(averaged over some zones and belts) , together with spatially resolved 
data at X = 8.11 ym and 8.45 ym (Orton 1975b) , indicate the presence of 
some obscuring matter at altitudes above the T = IbO^K level. The 
absorption spectrum of solid NH^ crystals fits this data quite well, with 
more absorbing material over zones than over belts suggested by the 
data. Nevertheless, a small amount of absorption seems to be required 
over the belts as well, mth an extent of only a few kms at levels with 
T ^ 145®K to 150"K indicated for this thin haze. 

One question of primary concern is the nature and location of the 
red chromophores which are responsible for the absorption of visible 
sunlight over the belts (or at least for the preferential absorption in 
the blue) . Upper limits to the location altitude can be inferred from 
the common observation that the red chromophores are sometimes overlain 
by time-variable white clouds. The relative heights ate apparent because 
the interface between white and red material is frequently convex outward 
from the white clouds. As discussed, the highest lying clouds which are 
possible on Jupiter are ammonia cirrus clouds at levels with T 2i 145 to 
loO“K. The red chromophores must therefore reside mainly at deeper 
levels '^^lth T > 145“ and P > 0.6 bar. Emission at X ^ 5y reveals infrared 
hotspots in only some regions of the belts (but in none of the zones) 
where the radi emerges from deep levels with temperatures up to 
T 300“K. The red coloration is anticorrelated with these infrared 
hotspots (Keay et al., 1973; Wes tphal et al. , 1974), which implies either that (a) 
the red chromophores reside above the level of some intermittent clouds 
which are opaque to 5y radiation and which reside somewhere above the 
T 2i 300“K level; or (b) that the red chromophores are themselves intermittent 
and opaque at 5y. The infrared hotspots are associated with visually dark 
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patches of bluish hue, probably connected with Rayleigh scattering rather 

I than with blue chroraophores, placing the blue regions at pressures of 

i 

! < 

j several bars, and T 'v 280®K (see Sagan 1971). To summarize: The red 

I chroraophores reside mainly in the troposphere (T > 145*K) > extending 

downwards from the tropopause (T 145'’K) to at most the T 300“K level. 

The condensation temperature of solid only slightly higher 

than the temperature (T 2l 145 °K) expected for the tropopause, so that 
the solid clouds reside in the upper layers of the troposphere. The 
pattern of patches of white clouds overlying the red chroraophores changes 
on a time-scale of a few weeks. This is comparable with the total tropo- 
spheric eddy diffusion time t^ in eq. (3) if 5 'u 1, but is also compatible 

with the longer total circulation time t , in Gierasch’s I 

® cxr j 

model, since time-scales for patches of horizontal extend d « D , (the 

P . i 

separation between zones and belts) are only of order (‘^p/® 2 b^*"cir* The postu- ‘ 

lated intermittent 5y clouds xmderlying the red chroraophores could in principle 
reside anywhere between the T ^ 145**K arid T ^ SOO'^K levels, but on 
I Gierasch’s model these levels in belt-regions represent a continuous : 

i downdraft of warming material which could not form NH^SH or H 2 O crystals* j 

Solid clouds can be present in belts to varying degrees ^ because the | 

I upper horizontal circulation path from zones to belts can include layers j 

; • ^ ■ ' • • ' - ■ I 

i just below the tropopause which carry crystals. On this picture j 

I the red coloration would have to be provided by red chromophores which | 

more or less coexist with clouds (P ^ 0.5 to 0.7 bar), but with i \ 

\ ^ " ■ j • 

j varying amounts of (i) high-level clouds with small NH^ crystals, (ii) red 

j. • - ■ 

I chroraophores, and (ill) clouds at loxrer levels with larger NH^ crystals. 

(c) Radiative transfer excursion i 

' Besides the red chromophores we are also interested in inferred . 

near UV-ehromophores: The measured albedo of Jupiter (Wallace et al. , 1972; : 
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Savage and Caldwell 1974) to solar photons in the near DV (X 'Vi 2000 A to 
3500 say) is quite low, with a minimum albedo 0.25 near X 2; 2800 X 
(after rising slowly towards shorter X, the albedo again seems to be lov 7 er 
for X < 2300 A) . Absorption by H 2 S , residing in the NH^SH clouds at 
lower levels, has also been considered (Prinn 1970; Lewis and Prinn.1970; 
Sagan and Khare 1971b; Lewis 1976). A computational problem is often 
presented by absorbers whose abundance decreases very rapidly with increasing 
altitude (because of condensation) ; we briefly review the relevant radia- 
tive transfer theory: 

Consider plane parallel geometry and assume the scattering is isotropic. 
Let X be the total extinction optical depth at some level (scattering plus 
absorption) , £5^ the single scattering albedo (ratio of scattering to extinc- 
tion cross-section) at the same level, and K = (3(1-S^)]*. If K is indepen- 
dent of X, the Eddington approximation is known even for slabs of finite 
thickness (Chandrasekhar 1960, Irvine 1975) and is discussed in our Appendix. 
The mean intensity J(x) as a ftmction of depth is particularly simple for 
a semi-infinite slab of constant ic (see eqs. A3 to A5 and eq. A17) . Unfor- 
tunately, K can be a very steep function of x as in the example in Fig. 2, 
which corresponds very roughly to absorption by H^S and Rayleigh scattering 
by H 2 of 2500 A wavelength photons in the work by Prinn (1970). Prinn used 
eq. (A3) for J(x) in the integrand of an integral which evaluates the 
absorbed energy layer by layer. Unfortunately, this use of eq. (A3) can 
lead to gross errors when K is varying rapidly, and a more accurate numerical 
evaluation is given in our Appendix for the example in Fig. 2. Fortunately, 
the correct result for the fraction f of the Incident energy which is ab- 
sorbed is fairly close to the following simple approximation, f^: Let x^ 

be the total optical depth at the level where K equals some predetermined 
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number, say K = 0.5 (in our examjile ^ 7.5). Let be the fractional 
energy which woiild be absorbed according to the Eddington approximation if 
one had no absorption at all for T < and a perfect absorber at T^. For 
diffuse incident illumination (instead of averaging over different angles 
of incidence) one has f^ == tO»75 x^ + 1] ^ if x^ » 1 (see eq, A16) . For 
the numerical example in the Appendix the value of f ^ corresponding to 
the choice ic = 0.5 for x^ fits the correct f quite well. 

(d) Hear PV i 

Solar photons in the near UV (X ^ 2000 A to 3000 A, say) are much 
more abundant than Ly a - photons , It is of interest to speculate about 
materials which absorb in the near UV to see (i) if their absorption can 
explain the low observed albedo in the near OV; and/or (ii) if the end- 
products of reactions initiated by such absorption can provide the red 
chromophores for the troposphere; and (iii) how much complex organic 
material can be so produced . 

The absorption of solar radiation with X 2300 to 2700 A by H 2 S 
in the upper layers of the troposphere has been discussed by the authors 
listed above (for X > 2700 X the absorption by H^S is too weak and photons 
with X < 2300 A are already absorbed by at higher altitudes). The 

1 O 

total flux in this wavelength range (diurnal mean) is about 1 x 10 photons 
—2 —1 

cm s and some fraction f of these photons are absorbed by H„S. We 

s ^ 

first present estimates if an ad hoe value of f^ = 0.1 is assumed [and ^ =1 
in eq. (3) ] . 
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Laboratory simulation experiments by near-ultraviolet irradiation 
of CH^, C^Hg, and H^S mixtures at 'v* 1 bar pressures have been per- 

formed, but without initial excess hydrogen fSagan 1968; Sagan and Khare 1971 a,b; 
Kiiare and Sagan 1973, 1975). Hydrogen sulfide is the primary photon acceptor 

in these experiments. Amino acids were produced in these experiments with 

-4 -2 

a quantum yield of 'v« 10 for molecular weights 100 ('^ 10 amu/photon) . 

The quantum yield of all organics in such experiments is of the order of 

100 times larger (of the order of 1 amu/photon) . Some red chromophores 

were abundantly produced* consisting mainly of polymeric sulfur, but including 

a substantial fraction of organic compounds. The overall mass absorption 

3 —2 

coefficient for blue light of these chromophores is h 'v* 10 gm cm . If 

such particles are produced in Jupiter at the rate of one 

.12 - 2-1 -11 - 2-1 

S-bond per absorbed photon (R 10 S cm s 'V' 5 x 10 gm cm s )» their 

column density,!?, over a few scale heights down to the thick water clouds, 

2 -45-2 

is needed. According to eqns. (3) and (8) , N 2RH 10 ' gm cm 

and the optical depth to absorption in the blue is t = hi? 0.03. This 

is almost two orders of magnitude smaller than the observed optical depth. 

Moreover the reaction of 8 sulfhydryl radicals to form Sg is poisoned by 

competing reactions, and the Sg quantum yield must be significantly < 1. 

In .'the laboratory experiments cited, organic molecules, but not 

polymeric sulfur, were produced from superthermal hydrogen atoms. In the 

presence of n LhouMJUul tlinuH mure Itj, (by muiibei )* Ibc qurtiiLiitii yield tiilgliL )jp pk.- 

pected to decrease by about the same factor, if energy loss of the fast 

_3 

atoms is mainly by elastic collisions, to 'v lO amu/photon. The total 
production rate would then be quite small, R 't* 5 o< ICT^^gm (which 

is still ten times larger than our estimated R for the mesosphere). 

org 

The corresponding optical depth fot organic chromophores is, then, many 

. ■ • ' ■ - . - ■ ■ - , . . i 

■■■■■■■ . . . . ••• i 



orders of magnitude below the observed value. In the laboratory experiments 
cited, large quantities of molecular hydrogen are produced during the ir- 
radiation because the end products are unsaturated and it may be that the 
correction factor from laboratory to Jovian situation is closer to 10 

_3 

than to 10 . It is also possible that somewhat longer wavelength ultra- 

violet light is atm effective in these syntheses, perhaps as long as 2900 A. 
Recent experiments <' (Ferris and Chen, 1975) have shown that 1849 A 

irradiation of mixtures of CH^ and la a 10 tl excess of H 2 produce a 
variety of organic molecules in very high yields, presumably by the genera- 
tion of hot hydrogen atoms from ammonia photodissociation. These results 
suggest that the effects of H 2 dilution on poisoning organic photochemistry 
may be slower than linear. Hence the contribution of organic polymers to 
near-UV photoproduced optical frequency absorbers on Jupiter may be sig- 
nificant. Polymeric sulfur is not pi educed by hot hydrogen atoms in such 
experiments and so the scaling from laboratory to Jupiter is more direct in 
this case. However, even h4rd, and even under the excessively optimistic 
assumption that every H 2 S photodissociation event leads to polymeric S, the 
implied optical depth falls short hy two orders of magnitude from matching 
the observed values; moreover, pure polymeric sulfur fits the observed 
optical properties of the Jovian red chromophores only poorly (Rages and Sagan 
in preparation). 

We have seen that a slow, deep circulation pattern is likely to extend 
all the way from the tropopause down to the level with T 300®K (Gierasch 
1976) with 5 'V 0.01 in eq. (3). Abundances of reaction products are propor- 
tional to and also to (f^/0.1) which was claimed to be as large as 4 by 
Prinn (1970). One might then expect H 2 S photolysis to dominate chromophore 
production, but this still seems unlikely for three reasons: (1) As discussed 
In Section IV (c) and in , iiii Appendix, Prinn* s (1970) estimate of f 'Vi 0.4 is 

S. 
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likely to be an overestimate since f = [0.75 T +1] is smaller. Thus 

o o 

f^^O.2 is more likely even in the absence of any scatterers (other than 

H„) or Other absorbers above the H..S levels. The smallness of f is also 
/ z s 

su gg ested by the fact that the observed albedo of Jupiter does not increase 
from X 2800 A Cohere H 2 S is Ineffective) to X 'v 2600 A (where H^S absorp- 
tion should already be strong) . (2) Since H 2 S cannot itself be the cause 
for the low albedo in this wavelength range, whatever absorbers are the 

cause must lower f further. (3) Most important, the systematic deep 
s 

circulation pattern mentioned before depresses H 2 S photolysis both in zones 
and belts; in the updraft, represented by zones, H 2 S is abundant in the 
NH^SH clouds but clouds are above these layers and prevent most of the 
solar radiation from penetrating down. In the downdraft, represented by 
belts, the clouds are absent (or at least diminished) but the material 
moving downward is practically devoid of H 2 O and H 2 S (w^hich condensed out 
on the previous upward journey). H 2 O and H 2 S are again abundant in deep 

f 

levels with T > 300“K but little of the near UV penetrates this deep. 

WHg absorbs solar photons of wavelength X < 2300 & which are 

O Cl 

somewhat less abundant than photons between 2300 A and 2700 A, On the 
other hand, should not suffer the shielding difficulties experienced by 
H 2 S because it resides in the highest cloud layer which extends up to the 
tropopause. NH^ should therefore be abundant in both zones and belts 
and even flows into the lower mesosphere at a non-negligible rate. It is 
therefore likely that most of the photons with X ^ 2300 & are absorbed by 
ammonia (Strobel 1973b), which is compatible with the observed decline in 
albedo from X > 2300 ^ to X < 2300 A. Hydrazine particles are also pro- 
duced by ammonia photolysis and may provide a suitable explanation for 
the low albedo near X 2800 & (Prinn 1974). Moreover, organic molecules 
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produced by photodissociation at X < 2300 K could conceivably be 

the principal source of red chromophores 3 if the quantum yields of such 
3 

reactions with 10 2 I dilutions are as high as the experiments of 

Ferris and Chen (1975) suggest. 

(e) Optical frequency photoproduction of chromophores 

We here considered two possibilities t If the circulation in the upper 

part of the troposphere is as slow as Gierasch (1976) suggests, then ^ < 0.01 

in eq. (3) and column densities obtained from eq. C8) are increased by two 

orders of magnitude; thus polymeric sulfur or complex organics might provide 

the red chromophores. But, as we have seen, there are serious problems with 

the contention that these molecules are produced by H^S photolysis. These 

problems are avoided if shorter wavelength photolysis of is an adequate 

source of complex organic chromophores, but this Is a questionable proposition 

considering the lower photon flux at X < 2300 1 , and whatever is the correct 

3 

value of the poisoning of organic photochemical synthesis in a 10 :1 dilution. 
If, on the other hand, C turns out to be of order unity, then ultraviolet 
photochemistry cannot in any case provide the observed chromophores at a suf-^ 
ficient rate. Thus, in both cases we are led to the possibility that 
photochemistry at optical frequencies, where more total energy is 
available, is involved in chromophore production. But no chemical 
bonds of even moderately abundant molecules can be broken by photons of wave- 
lengths longer than 3500 A on Jupiter, from which it follows that optical 
frequency photochemistry must be at least a two photon process. However, such pro- 
cesses appear to be exclusively biological. We therefore find ourselved led 
unexpectedly to the hypothesis that the Jovian chromophores are biological 
in origin and that there is an abundant biota in the Jovian clouds. 
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Let us assume that the optical chromophores are contained primarily 

vd.thin Jovian organisms in the upper troposphere — organisms driven by 

visible light photosynthesis, perhaps primarily in the blue, and utilizing 

the abundant methane, water and ammonia. Terrestrial algae (using CO^ 

instead of CH, ) under laboratory conditions (Lehninger, 1971) have a 

quantum yield for synthesizing glucose of about 4 amu per photon (in the 

visible) . The equivalent maximum production rate of biological material 

—9 —2 —1 

in the Jovian troposphere would he about R *= 3 x 10 gm cm s . With a mass 

3 -2 

absorption coefficient of k 10 cm , such organisms would yield an optical 

depth of T 2, slightly more than required by the observations. If the 

~2 

parameter ? turns out to be of order 10 , and the red chromophores are due 

to Jovian photcsynthetic organisms, R and/or K would have to be smaller. 

Thus Jovian organisms having metabolic and photosj'nthetic parameters 
typical of terrestrial algae are able to account for the optical chromophores 
on Jupiter. Such organisms must however be adapted to the Jovian environ- 
ment, and in particular must maintain a steady state population in the face 
of convective pyrolysis. Likely ecological niches for such organisms, as 

determined by the hydrodynamics on Jupiter, are discussed in the following 
sections* The following discussion can, however, be treated independently 

of the hypothesis that Jovian chromophores are biogenic. 


V. SINKERS AND FLOAIERS 


If there is a Jovian biology, we would expect it to fill a rich 
variety of ecological niches. The best terrestrial analogy seems to be the 
surface of the sea. Oceanic phytoplankton inhabit a euphotic zone near the 
ocean surface where photosynthesis is possible. They are slightly denser 
than seawater and passively sink out of the euphotic zone and die. But 
such organisms reproduce as they sink, return some daughter cells to the 
euphotic zone through turbulent mixing, and In this way maintain a steady 
state population. (The early stages of sinking move the phytoplankton 
from a region of depleted resources into a region of abundant nutrients, 
thus stimulating replication) . 

A more elaborate adaptation is provided by fish and other organisms 
with float bladders which use metabolic energy to maintain a habitat at 
suitable pressure levels. They are generaily not photosynthetic autotrophs, 
but rather heterotropfas living off organic molecules produced by autotrophs. 
A third ecological niche is filled by marine predators, one step fiirther up 
the food chain, which hunt heterotrophs . 

In the following discussion we will consider three comparable 
ecological niches on Jupiter: The primary photosynthetic autotrophs, which 
must replicate before they are pyrolyzed, will be described as sinkers. 

A second category of larger organism which may be either autotrophs or 
heterotrophs but which actively maintain their pressure level will be 
described as floaters. A third category of organisms actively seek out 
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other organisms; we call these hunters, although, as we shall sec, the 
distinction between hunting and mating under these conditions is not sharp. 
Finally, there is a category of organisms which live almost at pyrolytic 
depth. They are scavengers, metabolizing the products of thermal degrada- 
tion of other organisms. For our. purposes, these pyrolytic scavengers are 
identical with floaters. In the following two sections we estimate 
the growth time t required for an organism to double its 
mass. 

We first consider passive organisms with positive excess density Ap 

above the gas density p . Such organisms falling under gravity through the 

6 

atmosphere rapidly acquire a vertical drift velocity, w^, which depends on 
location and the size of the organism and increases with increasing Ap . 

There is a premium on keeping Ap^ and hence w^ small and we consider only 
organisms in the shape of thin, gas-filled balloons. Such organisms have 
been proposed briefly before both for Jupiter CShklovskii and Sagan, 

1966) and for Venus (Morowitz and Sagan, 1967). 

For simplicity we treat a spherical shell of outer radius a and 
skin- thickness d « a. In general we consider organisms filled With gas 
of the same temperature and mean molecular weight (2.2) as the surrounding 
ambient atmosphere. The lack of buoyancy is then controlled by the excess 
averaged density due to the skin, Ap = (3d/a)p where we assume the 

^3 

biological skin material to have unit density (l.gm cm ). 

For sufficiently small radius a of an organism, the Stokes formula 
in eqtt. (ll)holds (with ■+ 3d/a) for the drift velocity w^ and the drag 
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coefficient Cjj ^ 12v/aw^ is large. Let be a critical radius for which 
this relation gives Cjj = 1 (Reynolds number = 24) , so that 


d 3 vp„ 


®Re 


(17) 




For large organisms with a » (large Reynolds number) the flow around 

2 2 

the body becomes turbulent and the drag coefficient Cj^ (= drag force/0. Sira p^w ) 
is close to 0.5 (see, e.g. , Batchelor, 19 70 ) . In this regime 


■»d ^ 


(18) 


We shall have to consider separately organisms living below and above 
the tropopause. The troposphere provides conditions closest to terrestrial 
ones (similar temperatures and abundant water in liquid or vapor form) . We 
assume that the organisms can survive there over a range L of vertical height 
(L 5H, say), centered approximately on the level with T 'u 300°R. An impor- 
tant timescale is then a nominal drift-time t^ = l*/w^. We consider the skin- 

-4 

thickness d as an unknown parameter with d 10 cm possibly a typical value 

-5 ^ 

(biological unit membranes on Earth have d < 10 cm, but mechanical strength 
is required for the skin). Using eqns. (1) and (11) and evaluating 
coefficients at the level with T = 300®K, we find 


5H 


. T 1 cm 10~^cm , „ , „4 \ 

30D°K a d ^ 


for a < (V 
Re 'V 


10 ^cm 300 K 


h 


) ( 19 ) 


0.039 cm 


1 
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For B. > oil other hand, eqn. (18) gives a drift-i'ime which depends 

on d but not directly on radius a, 


^ ■ ^.,1 

d 5H 


T 


300“K 


,-4 


10 cm) 


2.5 X 10 sec 


( 20 ) 


The eddy diffusion coefficient in the convective troposphere is 


fairly large and we also have to consxder the contribution of eddy diffusion 
to the organism’s downward motion over a vertical distance L. We are 


Interested in cases where I. is larger than H and the effective circulation 


time t . for downward diffusion is then not given by L /K ; With L » H 
cir c 


and a slowly-varying function of height, the first derivative term in 


eqn. (7) can be neglected and eddy diffusion has the same effect as increasing 


w, to (w, + h”^) . The relevant time scale is then linear in L and is ap- 
d a 


proximated by the smaller of t^ and the quantity 


L LH „ L 

f- ■» _ . 

cir H c K 5H 
c 


5 

3 


300°K 


X 3.7 X 10 sec , 


( 11 ) 


where we have used eqn. (3) and assumed 5 'u 1, If d 10~^ cm, then eddy 
diffusion competes with gravitational fall only for microorganisms with 


a 30 )jm. Thus, Jovian organisms the size of small terrestrial protozoa 


and prokaryotes have typical times for falling through the troposphere to 
pyrolytic depths of one to two months. To maintain a steady state population, 
they must only replicate in that time scale. 
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In the Joirian mestjsphere the temperature is low (we assume T 150°K) 
and water is absent, but some biologies have nevertheless been envisioned 
for such environments (see, e^g., Pimentel et al«, 1966j Sagan, 1970). We 
shall see that the lower, denser layers of the mesosphere are the most 
advantageous for the growth of such organisms; eddy diffusion is negligibly 
slow here compared with the downward drift velocity w^. For organisms 
operating between a lower pressure P and the tropopause (pressure 1 bar), 
the drift-time L/w, becomes 


1 bar 1 cm 10 cm „ c 

— - — ^ X 9JU sec for 

F J a a 

> 

a < . ^ (io^ 0.044 m 

'v ^ d P j ^ 


( 22 ) 


For larger organisms, a > a^^ , the expression for large Reynolds number 
gives 


t 'V 


£n 


1 bar! 


10 ^cm 


1 bar 


1,6 X 10 sec (23) 


To obtain the pyrolysis time scale this number must be added to the 
appropriate tropospheric time scale given by eqs. (19) , (20) and (21). 

The expressions given above for t^ hold only for balloon organisms 
with ambient atmosphere Inside, For organisms capable of pumping gas, 
buoyancy can be achieved while maintaining pressure equilibrium by keeping 
the interior gas pure hydrogen with molecular weight p = 2.0 instead of the 
ambient hydrogen-helium mixture with ]i = 2.2, A pumped organism can thus 
float if Ap^ = (3d/a)p_ < O.lp . This requires radii a larger than a 

S O 
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threshold-radius for floating, given by 


^fl'^ 


^ 1 

T ^ 


^1 bar 

i_-4 

i50°K 



\ 10 cm 1 

1 


[ / 


3c 16 cm 


(24) 


Instead of maintaining a hydrogen interior, an organism could use a fraction 
of its metabolic energy release to heat its interior gas. The temperature 
difference AT between interior and exterior is determined by the metabolic 
rate and the thermal diffusion coefficient of the gas. Except at very high 
pressures, the thermal diffusion is fast enough to keep AT/T « 0.1, so that 

smaller values of Ap^ and larger radii than a^j^ are required. Buoyancy from 
a hot interior is thus likely to be less advantageous than pumping out helium. 
Note that Jovian floaters are macroscopic organisms (eq. (24)). Specialized 
floaters, with a range of organ systems, might have an effective value of 
d 1 cm as terrestrial reptiles and mammals do; in that case, floaters would 
have kilometer dimensions, a size within the resolution capability of the 
Mariner Jupiter/Saturn flyby imaging system. The existence of very large 
floaters is, however, limited by metabolic contraints described 


in the next section 


VI. GROWTH OF PASSIVE SINKERS AND FLOATERS 


We defer for the moment further questions of buoyancy and of powered loco- 
motion, and consider a range of metabolic niches. Organisms — sinkers or 
floaters — may be photosynthetic autotrophs, obtaining their free energy 
directly from sunlight. They may also be heterotrophs of which we distinguish 
two classes: (i) those which feed on organic molecules of non-biolog ical 
origin (e.g, , ethane, an ultraviolet photoproduct), which reach the organism 
by molecular diffusion; and (ii) organic matter in smaller organisms or their 
fragments which reach the larger organism because of the difference in 
vertical drift velocity w^ between predator and prey. Active hunters 
described in the following section are a related subcase. In cases 

we must estimate the growth- time t„ = [2 d£n a/dt]~^ and compare it with the 

drift-time t^ (or with t^^^ if that should be smaller). We shall see that the 

ratio t It, increases with the radius a of the organism; we are interested in 
gr d 

the maximum radius for which this ratio is unity. If is appreciably 

larger than the minumum size a . for a particular life-style, then a biology can 

rain 

be maintained by an organism fragmenting into many smaller organisms or dispersules 

before it reaches size a , i.e., before it has drifted do^wiwards to a pyrolytic 
— max 

death. If t^ » foi^ the dispersules, they can he circulated upwards, grow to 

Siax* complete the life-cycle. 

We saw in Section IV that the most favorable region for direct photosynthesis 
is probably the upper troposphere, just above the water clouds (with T 'V» 300°K, 

.. 40 , 
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P 'V 5 bar). Methane, atumonia and water building blocks are abundantly available 

there to a photosynthesizing organism and the growth-rate is controlled by the 

deposition of free energy from solar photons. Furthermore, the observed optical 

chromophores can be explained by such organisms at such a level. We saw that a 

_9 ^2 -1 

value for the production rate ofR‘'^3xlO gmcm s seems reasonable and we 
find for the growth- time 


= a/2 

gr - da/dt 


4d _ d 
R - 10-^ cm 


3 X 1 0 
R 


-9 


X 1.3 X 10^ sec 


(25) 


Equating this time with t^ in eq. (19) then gives for the maximum radius 

to which an organism can grow 


a 

max 


( 10 "^ 


cm/d) ' 


0.1 cm 


(26) 


Values of a , - 10*“^ seem perfectly feasible, so that one organism can produce 

mxn 

more than 10^ dispersules for which t^ « 4 x 10® sec The steady-state 

population of such organisms would then be controlled only by competition for sun- 
light, and not by gravitational fallout to pyrolytic depths. 

Consider next the diffusion onto the skin of an organism of some organic 
molecule with molecular weight u and abundance in the surrounding atmosphere of Pq 
( in gms cm“^ of atmosphere) . Ethane is the lowest mass (^ = 30) organic molecule 
which carries free energy (relative to methane, the prevalent carbon-molecule at 
thermal equilibrium) . The diffusion coefficient ]) of ethane in the Jovian H-He 
mixture is about 0,4 times the viscosity V of this mixture (Strobe! , 1973) and for 
heavier molecules B scales reciprocally as the cross-section of the molecule. Eon 
approximately spherical organic molecules we adopt 

D 0.4V (li/30)"2/3^ (27) 
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with V given by eq • (11) • If V > D ^ (aw^/12), the diffusion is almost indepen- 
dent of the motion of the sphere and the steady-state rate for diffusion onto the 
surface is Aira^ (p^ D/ a) . Let ^ ha the ratio of increased body-mass of the 
organism to the mass of ingested metabolites. Assuming unit density for the 
skin-material the growth-time is then 

” ad/2e pQ D. (28) 

At a mesospheric level with pressure P (in bars) , for molecules produced 
at a mass rate of R = at some higher level, eqs . (6) and (9) give 

Pq (]i$H/100) P S-nt The growth-time in the mesosphere is then 


t 'V 

gr 


a d 

1 cm id“^cm 


(u/30)^^^ 
£n (1 bar/p) 


r 5 X 10~^^- | 


X 1.8 X 10^ sec. 


(29) 


For ethane an efficiency factor e ^0.03 seems reasonable and 

R 'V 5 X 10“^‘^gm cm s~^ has already been discussed; for UV photoproduction of more 

complex organic molecules we may leave e 0.3 and R 5 x 10 gm cm s , so that 

2/3 

Re is the same hut the factor U makes heavier molecules less favorable. Replacing 
the slowly-varying factor (1/P) by 3 and Re by 1.5 x 10 we find from 
eqs. (22) and (29), 

a (IG”^ cm/d) (p/30)”^^^ x 0.07 cm (30) 

If is particularly large and/or d particularly small, an organism can grow 
beyond size a^^^ and the diffusion rate now depends on the drift-speed w^, given by 
eq* (18). We omit the complicated transition region a and make only order 

of magnitude estimates : There is a boundary layer in the flow around the sphere of 

thickness b a (4v/a « a, such that the flow-speed a normal distance y 
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from the body Is 'v» (y/2b) w^. With D < v the typical stand-off distance for dif- 
fusion onto the body-surface is bjj ''' 2b (D/v)^^^ and the diffusion rate is 
^ 4ir a^ (pQ D/b])), so that 

tg^ ^ (2d/ePp) (a/Dw^)^^^ Cv/D)^^^. 

For the mesosphere this gives 


gr '^1 


1/2 


cm 


,_d ^,3/4 ifear IM 

'l0“4 cm' ^ P ' in (1/P) 


The same basic formulae apply in the troposphere, but p^ is now much smaller 
for the same value of $ because of the larger eddy diffusion coefficient, K^. On 
the other hand, we have seen that photosynthesis may be going bn there and this 


may release complex organic molecules at a faster rate 0. Substituting p^ = RH/K^ 

and using eqv (3), we find (for a < a^g) 

^ a d / T Uv2/3 r 3 x 10'“^ 0.3i „ .rJ 

1 cm io~^" cm ^300^K 30^ ^ $ e ^ ® ^ 


Comparison with eq. (19) shows a ^^y to be of order (10*"^ cm/d) (30/y)^^^ 0.013 cm, 
which is just macroscopic. 

If small tropospheric photosynthetic organisms have radii a ~ 0.01 cm, their 

s 

drift speeds are even smaller than the eddy icirculation speed and their abundance 
is still given by p^ - RH/Kg. A larger organism of radius a, drifting downwards 
with speed w^, can grow by coalescing with such small organisms if the encounter 
satisfies certain conditions. These conditions are similar to those governing the 
coalescence of raindrops (Mason, 1972; Rossow, 1975). In analogy with eq. (13) 
one finds, independent of a_, and provided d « a^ < a, the requirement that 
a > (10“^ cm/d) (T/3Q0°K) 0.002 cm. If these conditions are satisfied, the larger 
organisms incorporate mass from the smaller ones at a rate of Att Pq that 


(31) 


(32) 


(33) 
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A/ (4d/R) (K^/Hw^j). If a ^ a, 0.04 cm, then 


t ^ (.. ..4 ,, . [3 ^ ip:^. %3] 2.5 * in'* sec. 

gr '‘10 ^ cra*^ '300^ # ^ 


Neither tg^r ao^ t^j in eq. (20) depend explicitly on radius a. An organism can 
then grow indefinitely as long as (d/10~^ cm) $ (T/300*^K)^^^. 

Thus, provided the membrane thickness stays smaller than about 1 ym, such 
passive sinkers can grow by coalescence with smaller organisms to sizes larger 
than a^^ as given by eq. (24). This provides one evolutionary pathway from 
sinkers to floaters. However, the condition on membrane thickness is very restric- 
tive. Floaters of this sort are unlikely to be much more stable than soap bubbles 
and the organ systems which large organisms require could not be accommodated with 
the thin skin thickness. However, as we shall see in the following section, the 
hunting lifestyle provides a way out of this dilemma. 

To summarize the situation so far: Sinkers can grow and produce reproductive 

dispersules in a stable steady state life cycle if they either (a) are mesospheric 
heterotrophs or photoautotrophs, or (b) are tropospheric photoautotrophs. Passive 
floaters exist only if (1) their skin can be made particularly thin, and (2) a 
significant fraction of the smaller photoautotrophs is available as food. 



VII. POW]EEED LOCOMOTION AND COALESCENCE 

We saw that the efficiency of utilization of metabolites, diffusing onto the 
surface of a freely falling organism, decreases with increasing mass of the molecule. 
However, for even more massive metabolites and organisms approaching each other with 
some relative drift velocity, growth by coalescence is again moderately efficient. 

The groxi?th-rate by coalescence can be increased greatly if organisms are capable of 
powered locomotion and can steer toward each other. We shall not discuss the 
predator-prey relationships of a strict hunting doctrine, but rather only an ecology 
in which an organism approaches another, coalesces into one larger organism 
which approaches other organisms, etc. It is clear that this hierarchical mating 
doctrine car also serve the traditional terrestrial function of mating, namely, the 
exchange of genetic material. A strict hunting doctrine in which one of the partners 
in such a coalescence is non-cooperative will be less efficient than the doctrine 
discussed here. For convenience, however, we will describe both doctrines as hunting. 

In the following discussion we will assume that only one partner in such a mating 
event is active in the search; for this reason mating, hunting prey, and hunting 
abiological organic molecules are not extremely different. The doctrine in which 
both organisms are active in seeking partners would be more efficient than the ones 
described below. 

For simplicity we assume the presence with number-density n, of organisms with 
radii of order Sl-^ at some upper inj ection— level z = The organisms coalesce and 

grow as they drift downwards, so that the size distribution changes With height z. 

Let n (ai, z) and p(a; z) .= (4 f d) n (a, z) be their steady-state number-density 
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and Kiass“density, respectively, per logarithmic interval of radius a. In principle 
we have a relation like eq. (7) separately for each value of a, with the 
right-hand side representing the flin: from one size range to another, tfe shall make 
only order of magnitude estimates for a(z), the mean radius of organisms at level z, 
and for n(z) = / djilna n(a, z). 

We assume that the organisms are sufficiently large that w,(z) , the drift-velocity 

a 

for radius a (z) at height z, is much larger than K/H in eq. (7). Multiplying by 
mass (assuming as usual unit density for the skin-material) and integrating over sizes , 
we find 

E4ira~^(z) d] Wj(z) n(z) - e(z) R (35) 

In this relation R is the constant rate in gm cm~2 at which mass in the form of 
organisms enters the upper level z = Z-j^; e(z) is a slowly-varying factor (slightly 
smaller than unity) which takes account of any inefficiency in the utilization of 
biological material, due to waste or expenditure of free energy. We shall have to 

evaluate the growth-time t (z) S [2 d^n a/dt]~^ for organisms with a “v a(z) at 

. ■ 

level z. The function a(z) is then given in terms of tg^ by 

dJln a/dz - [2 w^ (z)] ^ (36) 

For the integrated column-density N(a) = fdz n(a, z) we also find N(a) 'u 2 t w^ n(z) 
where z is such that 3(z) = a. The total mass contained in organisms in different 
size ranges is then, given by 

(4ir a^d) N(a) 2 tg^(a) e(z) R. (37) 

Assume that an organism, while drifting doxraward with speed w^, can expend 
enough mechanical energy to acquire a horizontal velocity component of up to 9 
where 0 is some angle less than 1 radian. Its vertical speed relative to other 


L. 


i 
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similar organisms is about 1/2 w and in time t a volume (0^/6) (w, t)^ is accessible 

d u 

to a searching organism, When n(z) times this volume is unity, the organism has a 
good probability of coalescing with a mate and doubling in mass. With n(z) given by 
eq. (35), the growth-time is then 


t Oj 

gr JJ.n2 


L. ( 3TTa^ d s 
n2 eB.0^ w^^*^ 


,1/3 


(38) 


If a > we have, independent of d both for the troposphere and the mesosphere, 

(39) 


t [1 (_£_)^ 4Z5fK 

’^gr '^1 cm'^ T 


P 3 X 10"® 0.3,1/3 


For such organisms the troposphere is quite favorable even in the absence of photo- 
autotrophs as prey (R 5 x 10*”^^ g cm“^ s”^) ; for instance, with T 'u 300®K and 

d 'v, 10“^ cm, we have a_, 'V' 6 cm and a (6/0.002) 6 cm, so that buoyancy can be 

1 1 max 

-9 _9 — 1 

achieved with quite small values of 6, For R ''j 3 x 10 g cm”* s , 

a^j^ 'v> 6 cm requires only a^ ^^ (6/3 x 10“^) 6 cm; microradian maneuverability suffices. 

The maximum sizes of organisms embracing the hunter doctrine as a function of the 
maneuverability, 0, and the source of food is given in Table II. We see that very 
large hunters are permitted. 

To derive an explicit relation between typical radius a and vertical level z (a), 
we must Integrate eq . (36) with t given by eq , (39), The result depends on 

gT 

whether 9 is constant or a function of a. Neglecting the variations of 6, w^ and e 
gives (for a much larger than the original size a^^ at level Z^) 


- z (a) ■= 3 w^ (z) tg^ (a) 


(40) 


We still have to calculate the mechanical energy expended by an organism in 
maintaining its horizontal hunting speed 6 w^. In pure downward drift the rate 
of energy dissipation equals the rate of gravitational energy release mg w^. The 
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TABLE II 

Maximum Hunter Sizes 

Source of Metabolites 

Abiological organic Photoautotrophs 
matter 

e 

~3 

2 X 10 radians 6 cm 40 m 

10“^ 3 m 2 km 


1 


30 m 


20 km 
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equivalent rate fpr maintaining the horizontal speed (with 9 < 1) is at most 
0 mg (^<rith judicious sailing and soaring it could be appreciably less). The 
relevant quality factor is the ratio &E/Am of mechanical energy expended in find- 
ing a partner to the increase in mass of biological material. (Biochemical energy 
storage is of order 10^^ erg/gm and a mechanical quality factor 10® erg/gm is a 
very modest assumption.) The heat of combustion (full oxidation) of mojecular 
hydrogen is about ten. times that of glucose and more than three times that of such 
fats as stearic acid. It Is not the unavailability of reduced compounds but rather 
the unavailability of oxidized compounds on Jupiter which probably sets limits on 
the efficiency of internal metabolism of hypothetical Jovian organisms. We have 


M ■JC 0 g w , t , 
Am d gr 


(41) 


which depends quite weakly on temperature and pressure. For the upper troposphere 

d 3 ^ ^ ere (42) 


^ (I7i4— Ce 5 X 106 SS6. 

Am 1 cm 10 ^ cm R e gm 


represents quite a modest power expenditure except for the largest organisms. A 100 m 
radius hunter of photoautotrophs with 0 = 0,01 radian excursion capability and a 1 cm 
skin thickness requires a quality factor of 10^^ erg/gra. The same quality factor is 
required by a 1 m radius hunter of abiological organic matter with 0 = 0.03 radians 
and d = 1 cm. 

Even smaller power expenditures will be required if, as would be likely, evolution 
selects organisras with sensory systems which increase the efficiency of hunting food — 
such as, e.g., optical sensors for detecting chromophores . Likewise, acceleration 
sensors capable of detecting convective "thermals” would ease the restrictions on 
floater sizes and energy expenditures; and habitats near the boundaries 
between upward and downward convective regimes would improve the cost 
of Jovian biology. 

The path from sinkers to hunters to floaters is an evolutionary path. It pro- 
vides a possible sequence for the evolution of floaters which, as we have seen, would 
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be difficult to understand directly from sinkers. Once floaters are established, 
they may be able to replicate without passing through the sinker stage at any point 
in their life cycles. For this to occur, of course, the minimum radius of their 
dispersules must be larger than 

While we have distinguished passive sinkers, hunters, and floaters 
in this discussion, it is clear that a single highly evolved Jovian orgatiism 
might partake qualities of all of these lifestyles ““ for example-, at different 
stages in the life cycle. Indeed, it may be that only after "he hunting 
doctrine is adopted by sinkers can growth to sizes large enough for the 
evolution of floaters become possible. We have in this discussion made no 
distinction among various locales on Jupiter; but it is clear that some 
locales -r- th^ Great Red Spot, for example — may be more favored than others 
because of higher abundances of organic molecules , prevailing updrafts 
or other reasons. 

Among other objects in the outer solar system which possess atmos- 
pheres, only Titan and Saturn exhibit colorations similar to those on Jupiter. 
Because of the absence of convective pyrolysis on Titan, the chromophores 
there can quite readily be ultraviolet photoproduced organic molecules 
(Sagan, 1973). However, the physical environment of Saturn is very similar 
to that of Jupiter in the respects relevant to the arguments of this paper 
and we, therefore, tentatively postulate an airborne biota on Saturn as well. 

The test of these ideas is, of course, observational. Flyby, orbital, 

and entry probe spacecraft each can perform significant tests of these 

ideas. For example, entry gas chromatograph/mass spectrometers for the 

near future are anticipated to have detectivities of organic matter as 
. _io ' _3 

good as about 10 gm era « The predicted steady state density of complex 
organic molecules on Jupiter is and is roughly estimated, using the 

numerical values of this paper, both for the lower mesosphere and upper 
troposphere and for abiological and biological sources of organic matter in 
Table XII. While most biological organic matter will be in organisms and 


TABLE III 

Steady State Density of Organic Molecules 

Abiologieal 

Lower Mesospliere 2 x 10“^1 

’“1 s 

Upper Troposphere 2 x 10 


(g cni-3) 

Biological 
9 X 10“® 

9 X 10"^^ 
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not in free molecules, the densities predicted in Table III are suf- 
ficiently large that a high sensitivity entry mass spectrometer may 
be a significant test of the biological hypotheses of the present 
paper. Because of the slox -7 convecuive velocities of the mesosphere, 
such instruments should be capable of working at such pressure levels. 

The possible existence of indigenous Jovian organisms is also 
relevant to the question of sterilization of spacecraft intended for 
entry into the atmosphere of Jupiter. Even if the ambient environ- 
ment is inconsistent with the growth of microbial contaminants from 
Earth, the internal environment of Jovian organisms may be much more 
hospitable. The replication of terrestrial contaminants in the 
Jovian clouds also depends on the availability of trace elements. We 
have seen (Section III) that indigenous Na is probably missing from 
the upper troposophsre; and that indigenous Mg must certainly he 
missing. Magnesium ions are essential for nucleic acid replication 
and a wide range of other biological functions; if Mg is missing 
from the clouds, the likelihood of terrestrial biological contamina- 
tion of Jupiter becomes nil. However, exogenous sources of sodium, 
magnesium and other trace elements, particularly from micrometeorites, 
while small, canhot be neglected (Section III). Consequently, it 

seems judicious not to exclude prematurely the possibility of bio- 
logical contamination of Jupiter by terrestrial microorganisms. 

Acknowledgment 

This research was supported by NASA Grants NGR 33-010-101 and 
NGR 33-010-082 and NSF Grants MPS 72-05056-A02 and MPS 74-17838. 


53 


We are grateful to Drs^ R. Danielson, B. Draine, 
P. Gierasch, K. Hare, N. H. Horowitz, D. Hunten, 
J. Lederberg, J. Lewis, R. Prinn, W. Rossow and D 
interesting discussions. 


F. Dyson, M. Eigen, 

A. Ingersoll, W. Irvine 
. Stevenson for 


r 


■J 1 -I 


APPENDIX 


Consider a plane parallel slab of total optical thickness and 
let I (T,]i) be the intensity of radiation at depth T with |j the cosine 
of the angle between the propagation direction and the forward normal. 
Assume isotropic scattering and let be the single-scattering albedo, 
which may be a function of T, The equation of radiative transfer 
(Chandrasekhar 1960; Irvxne 1975) then reads 


U ^ I'(T,y) « 


dy'I* (T,y )+ S^e 
-1 


(Al) 


for an incident beam of unit intensity with direction cosine where I' 

— T /it 

is the intensity of the diffuse radiation and e ° that of the attenuated 
incident beam. The Eddington approximation can be obtained by substituting 
a function of a single variable for I’ (T,y) when 0 < y < 1 and another 

function I’ (x) when -1 < y < 0. Let 


1C = 13(1 - 


(A2) 


and let J’ (t) be the average over all y of I(T,y). 
semi-infinite medium (x <») one finds 


J'(T) 


-KX 


ae 




5 


For constant K and a 

\ 


l+(2/3)y“^ , , ^“o 

^ = ilT(2/3)K “ ^ ^2 


> 

(A3) 

/ 


The effective reflection coefficient R = 1 - f, where f is the fraction 
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of “the energy flux that is absorbed » is given by 

^o 

^ “ (1 + K\i^)[l + (2/3)lc] 


(A4) 


Generalizations for non-isotropic scattering are discussed by Chandrasekhar 
(1960), Prinn (1970) and Irvine (1975). For isotropic scattering Prinn 
used eq, (A3) for normal incidence = 1) in his numerical work. For 
general y^ but with K « 1 eq, (A4) reduces to 


f (y^) = 1 - R(y^) k( 2/3 + y^) 


(A5) 


For many purposes one is interested only in results averaged 
over a whole daylight-hemisphere; i.e. , averaged over y^ from 0 to 1. 

Such averages are represented by the single problem we shall consider 
here, where the incident intensity equals unity for all y between 

0 and 1. In the spirit of the Eddington approximation we replace the 
problem by one for intensity I(T,y) without explicitly isolating the 
attenuated incident beam: In analogy with eq. (Al) we have 


i(T,y) - 


-I(T,y) + jdy’I(T,y») 


-1 



(A6) 


which leads to the coupled equations 


- -j(T)[i - a^(T)i , 



(A7) 


/ 
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where 


J(T) 




I dyI(T,u), 
-1 


F(t) = h 


f duviI(T,u) 


-1 




> 

2 

K(t) = % j dyii I(t,u) 

J 

We again make the approximation (but here for all the radiation) 
I (t) for u > 0 and I(T,y) = I_(t) for y < 0, so that 


(A8) 


that I(T,y) 


J “ 3K = + I_) , F = %(I^ - I_) (A9) 

In this problem, any explicit consideration of the incident beam is 
simply replaced by the boundary condition that “ 0) = 1 and the 

effective reflection coefficient R = 1 - f is simply given by I_(t = 0). 

For general and with k(t) defined by eq. (A2) , the use of 

the Eddington approximation (A9) iii eq. (A7) leads to the eikonal equation 


dT 


i+Cir) = 0 


and to the further relation 


2k^I^ + 1_) = 3^ (I_ - I^) . 


(AlO) 


(AH) 


We first solve these equations for the case of constant K for a slab of 

total thickness T : We assume that the effective reflection coefficient 

o 

R S I (t ) /I, (t ) is given as one boundary condition, the other boundary 

O —0+0 

condition for the incident beam being I (0) - 1. The general solution of 
the coupled equations (AlO) and (All) Is 


\ 


I+(T) = a[(3 ± 2ic)e"'^'' - y(3 2K)e*^* 


(A12) 


with a and Y disposable constants. The requirement of I /I, “ R at the 

“T O 

second boundary (T = T^) leads to 


(3 - 2tc) - R^(3 + 2 k) 
(3 - 2k) - R^(3 - 2k) 


(A13) 


and the condition 1,(0) = 1 then gives a. The final expression for the 
intensity J (averaged over direction cosines) is then 


J(T) = - Y« . 


k(T - 2T„)j 


(A14) 


(3 + 2k) - (3 - 2k) Ye 


Another quantity of interest is the effective reflection coefficient R^ 
of the whole slab at the surface of incidence, which is I (0), 


^ « (3 - 2k) - y( 3 + 2K)e~^'^'^P 
^ (3 + 2 k) - y(3 - 2k) 


(A15) 


Two limiting cases are of interest for the solutions given by 
eqs. (A12) to (A15). Consider first the case with R^ = 0 in the limit of 
K ->■ 0 with kept finite. In this case eqs. (A14) and (A15) reduce to 


. 3(t - t) + 2 

' = 1 \ - 3T-T« ■ - 3° 

o o 


(A16) 


The other case is the limiting case of -»■ “ for any constant, non-zero 
K. Independently of the value of R^, this case gives 

4k 3e"‘^'^ 

f - 1 \ - 3 + 2k ’ ^ 3 + 2k 


58 


When < « 1 the value of f in eq. (A17) agrees with the average over 
0 < < 1 of eq, (A5) which is f = (4lc/3). 

No exact solutions of the eikonal equation (AlO) exist for general 
fc(T), but if tc varies sufficiently slowly, one is tempted to use the WKB 
approximation, for a semi-infinite medium, only one of the two WKB 
solutions survives and we have 

T 

J<t) oc expt-| dx*K(T')l (A18) 

o 

Note that substituting a varying k(t) into eq, (A17) (or into eq. (5) of ^ 
Prinn 1970, which is the equivalent of our eq. (A3) for normal incidence 
instead of diffuse illumination) is not a good approximation to eq. (A18) 
if K varies appreciablji' even if the variation is slow. Furthermore, the 

2 

requirement for the validity of the WKB approximation is that (dK/dt) « K , 

which is not satisfied by practical cases for k(t), such as the example 

in Fig. 2. An inelegant but reliable numerical method proceeds as follows: 

Subdivide the total slab into a large but finite number of thin slabs 

(with At 0.3, say) and replace (c(x) for each thin slab by its average 

value over this slab. Start with the last thin slab, assume R = (3 - 2 k)/ 

o 

(3 + 2k) and evaluate eqs. (A13) to (A15) with the slab- thickness x^ replaced 
by Ax. The value for obtained from eq. (A15) for the last slab is then 
used as R^ for the penultimate thin slab and J and R^ are again evaluated, 
and so on Iteratively until the first thin slab. The solid curve label!) ’ 

J is the result of such a numerical solution for J(x) and the effective 
reflection coefficient R^ for the whole medium (obtained at T = 0) is R^ jv 0.85 
The dashed curve is the result for J(x) in a different problem — no absorp- 
tion for X < T^ and complete absorption for x > x^, with x^ 7,5 (the value 
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of T for which K = 0.5 in the actual problem). Eq. (A16) gives the solution 
for J(t) for this problem and also yields 0.85. Eq. (A16) for should 

be a fairly good solution for any k(t), with the value where K 0.5j 
as long as » 1 and k is small for T appreciably less than T^. 

The accuracy of the numerical solution described above is limited 
(apart from computational errors) by errors inherent in the Eddington 
approximation, including the use of a constant coefficient in eq. (All). 

For a single slab this approximation is excellent when K « 1, although 
other two-stream approximations are better (Sagan and Pollack 1967; Irvine 
1975) when « 1. For the problems with » 1, discussed here, k is 
small near the incident surface and increases only slowly with T (at first) . 
The results for J(t) and R^ are most sensitive to layers with moderately 
small T and are insensitive to inaccuracies in the Eddington approximation 


at large x. 
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Fig. X: A schematic plot of temperature , T, and eddy diffusion 

coefficient, K, versus pressure, P for the Jovian 
atmosphere. The dashed portion of the curve for K 
corresponds to Gierasch's (1976) model for a deep 
circulation pattern driven by latent heat effects. 

Fig. 2: A model radiative transfer calculation for an assumed 

form of the parameter, tc, as a function of total 
optical depth, T. Numerical results for the mean 
intensity J are given in the solid curve. The 
dashed curve is an analytic approximation, discussed 


in the text. 


Ammon ia*ClPLidS“ 

Water'-CJouds ■ 


Ammonium Chloride^ 


Silica 



ibl 

Chromophoresr 
H2S Photolysis 


Pyrolysis 
of organic- 

compounds 
1 ' 


000 2000 5000 

) 









I 

OJ 

i 

0.5 

0.3 

0.2 

O.iS 

OJO 

0.07 

0.05 

0.03 

0.02 


